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Finding a consistent way to stabilize the various moduli fields which generically appear in string 
theory compactifications, is essential if string theory is to make contact with the physics we see 
around us. We present, in this paper, a mechanism to stabilize the dilaton within a framework that 
has already proven itself capable of stabilizing the volume and shape moduli of extra dimensions, 
namely string gas cosmology. Building on previous work, which uncovered the special role played by 
massless F-string modes in stabilizing extra dimensions once the dilaton has stabilized, we find that 
\Q • the string gas cosmology of such modes also offers a consistent mechanism to stabilize the dilaton 

itself, given the stabilization of the extra dimensions. We then generalize the model to include 
, D-string gases, and find that in the case of bosonic string theory, it is possible to simultaneously 

stabilize all the moduli we consider consistent with weak coupling. We find that our stabilization 
' mechanism is robust, phenomenologically consistent and evades certain difficulties which might 

, previously have seemed to generically plague moduli stabilization mechanisms, without the need for 

Pj^ . any fine tuning. 
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' The question of how to stabilize the multitude of moduli fields that arise in any attempt at obtaining 4-d effective 
t-H , physics from stringtheory, is at present, still an open one. Although there have been many notable attempts in the 
last few years (see [l| and references therein), many issues remain open. The majority of these attempts are crouched 
in the framework of low energy effective field theory, where either ad-hoc moduli potentials, potentials generated 
by higher order (e.g. in a') effects, or potentials generated by fluxes are used to stabilize the various moduli fields. 
In spite of the various successes encountered by these attempts, they come up against some quite general problems, 
i-Q ■ consistency with late time cosmology not least among them. A specific problem that accompanies moduli stabilization 
schemes that use potentials to fix the moduli fields, is the fact that any such potential will generate a cosmological 
constant at its minima, which introduces a host of other theoretical problems [2( . 
^ 1 String gas cosmology (SCG) 00 QUI, on the other hand, has to date focused on a narrow subset of the moduli 
problem- that of stabilizing extra dimensions. From its origins as a mechanism in which a 3+1 dimensional universe is 
generated by the dynamics of strings in the early universe |3J |J] , it has become clear that the effects of a string gas can 
stabilize any number of extra dimensions at the string scale in a way that yields a late time scenario that corresponds 
to our present day universe 001' provided that the dilaton has been fixed through some external mechanism. Key 
to this stabilization mechanism, and to the consistency of the SCG program in general, are the properties of a gas of 
massless string modes. It is the goal of this paper to report on even further-reaching consequences of massless string 
gas cosmology, which evidently contains a natural mechanism to stabilize all moduli fields for a particular toy model 
(bosonic string theory compactified on a torus). 

We will shortly provide an a priori motivation for focusing specifically on massless string modes if SGC hopes 
to be theoretically and phenomenologically consistent, after which we will begin in earnest our attempt at moduli 
stabilization using these modes. Before we do any of this however, we should offer a brief (and neccesarily incomplete) 
historical review of the SCG program, in order to give us perspective on why we feel it is rather strategically placed 
to solve the moduli problem. 

The Brandenberger-Vafa (BV) mechanism (proposed in 0), posits that the observed dimensionality of macroscopic 
spacetime arises from the dynamics of strings in the early universe. Such a mechanism was put on firmer footing 
in 0, where it was reasoned that dilaton gravity (i.e. a dynamical dilaton) was neccesary for such a mechanism to 
work. Such reasoning was based on the observation that any positive energy density and negative pressure (which 
accompany winding string states) in gereral cause the expansion of spacetime in general relativity, rather than its 
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contraction or stabilization, whereas for dilaton gravity, negative pressures generically cause contraction. In the years 
that followed, much work has gone in to putting the ideas of 0] on a firmer footing. Among the more notable results 
has been the demonstration that the BV mechanism still holds once we include branes of various dimensions as well 
as strings as the fundamental degrees of freedom ||. Since then, it has been shown that, in the context of dilaton 
gravity, string gases can effectively stabilize extra dimensions at the string scale |15| , with a resulting isotropy of the 
non-compact dimensions 0] . More recently, it has been shown that string gases can stabilize any number of extra 
dimensions even after the dilaton has stabilized @E3j m a wa y that yields a resulting late time FRW cosmology, 
and which has a consistent phenomenology (which we will uncover for oursleves further on). There is beginning to 
emerge, a surprisingly complete cosmology that results from having string gases pervade the universe, which we hope 
will eventually provide us with testable predictions on the nature of dark matter and on effects on the microwave 
background Uj, in addition to providing a framework in which one can potentially model non-singular bouncing 
cosmologies |l2l . Even more recently, it has been shown [T^ that, in addition to stabilizing the moduli corresponding 
to the radii of the various cycles of a toroidal compactification, the shape moduli are also stabilized. Emboldened 
by this, it might appear that if SGC can stabilize the dilaton itself, whilst preserving the stabilization of all of the 
moduli associated with the compactified dimensions, then it might serve as an example of how moduli stabilization 
in more complicated string theory settings might be effected. 

As mentioned previously, the workhorse for this paper is going to be a gas of massless string modes. Exactly 
why this should be, will need no further justification if we can demonstrate moduli stabilization using these modes, 
but attempting an a priori answer to this question provides a useful check on the theoretical and phenomenological 
consistency of what we are doing, and answers some of the criticisms that have previously been raised against SGC. 

A. The Case For Massless String Modes 

Consider first the idea that a gas of strings can stablize the various moduli fields that might appear in our particular 
string compactification. An immediate objection one might have is that there is no way to do this without severely 
overdosing the universe. That is, since all string states have masses which go like the Planck mass (if we equate this 
to the scale of stringy physics): 



m „ — ~ M pl N, (1) 
Va 

where N is a number of order unity for the lowest lying string states, and M p i = 2.43 x 10 18 GeV is the 4-d Planck 
mass, one might expect a fluid of such strings to be phcnomcnologically ruled out. However, if we consider a gas 
of such strings which are massless, then one has a shot of evading phenomenological bounds. In and ^} it was 
shown that a gas of massless string modes can effectively stabilize any number of extra dimensions at the string 
scale without overdosing the universe, whilst simultaneously providing a strong enough stablization mechanism that 
evades fifth-force constraints 29] . However, there is another important theoretical reason for us to focus exclusively 
on massless string modes. We will very shortly be considering dilaton gravity in the presence of string sources, and 
accordingly will be in the regime of a low energy effective theory described by the action: 

S = 2^2 / d D+1 x^—ge- 2 * (r + 4^*0"* - -^H^ X H^ . (2) 

Where $ is the dilaton, is the D + l dimensional metric tensor, H^x the field strength of the two form field B^ v , 
and Kq is the dimensioned normalization of this action. This is specifically the 'massless background' |23j[4| of the 
low energy effective fields, as it is only valid for curvature regimes such that: 

a'R c « 1, (3) 

where R c is the characteristic curvature of the background in question, and a' is related to the string tension as 
T = l/27ra'. As long as such a condition is satisfied, the curvature scale of the background is much bigger than the 
string scale, and massive string states will not be created. Conversely, if we wanted to introduce massive string states, 
then we would have to include higher order corrections in a' . Hence if we are to introduce a gas of strings as the 
matter content for the action above, then it should only consist of massless string states, as this action describes a 
consistent background for these states alone. 
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A potential concern that one might have however is that, naively, one might violate JSJ if we are considering 
compactifications of the extra dimensions that are at the string scale. However, it turns out that when we consider 
the spacetime dynamics of such a gas on a toroidally compactified background, which eventually results in the 
stabilization of the extra dimensions and the dilaton, and compute the components of the Ricci scalar throughout the 
evolution of the system, we will find that the condition © is satisfied to an excellent degree. Hence our framework 
is consistent throughout, and we will be justified in using the low energy effective action (i.e. dilaton gravity) as 
the framework for our physics. Further issues concerning aspects of the string spectrum which we will be using are 
considered in the appendix, where we derive the energy momentum tensor for a gas of strings the backgrounds of 
interest for this report. There, we find that a string gas contributes an energy momentum tensor that can be derived 
from the following action: 



S = - J d D+1 Xy/g^E\g a ,<l>], (4) 

where E is the internal energy of our string gas, which depends on the spatial part of the metric and the dilaton 
field. From this, we immediately infer that SGC has several compelling advantages over the effective field theory 
approach to moduli stabilisation. We see from the form of the above action, that one does not have the usual fine 
tuning problems associated with introducing potentials for the moduli fields. Such a potential would contribute to 
the energy momentum tensor from the action: 



S = - I d D+1 x^gV(4>), (5) 

where one would have to fine tune the minimum of the potential to vanish at the stabilization point to avoid an 
unacceptably large De-Sitter phase of evolution after stabilization. We see from Q that the coupling of the string 
gas to gravity is rather more complicated than just the usual yfg x V dependence, and the 'potentials' one derives 
from the moduli fields from do not require their minima to be fine tuned at all, and in fact, quite generically result 
in Freidman-Roberston- Walker (FRW) expansion of the non-compact dimensions^!] 30]. In this way, we evade the 
negative results of Giddings |2j concerning radial moduli stabilization. 

However from a purely phenomenological point of view, SGC has yet another advantage over the field theoretical 
approach to moduli stabilization: consider the fact that any potential one introduces to stabilize the moduli will, from 
the prespective of the 4-d effective theory, cause the modulus in question to have a mass given by: 

m 2 ~ V"{4>), (6) 
whereas the contribution to the energy density of the universe of such a stablization potential will go like: 



P~V(4>), (7) 

which, if we are close to the minimum of the potential, is also proportional to m 2 . When we try to make contact 
with present day observations, the simultaneous requirement that the mass of the moduli field be large enough so 
that we evade fifth-force constraints, whilst small enough so as not to overdose the universe becomes very difficult to 
satisfy, and introduces its own host of fine-tuning problems. The case for string gases is markedly different. We will 
see further on that the 4-d masses of the moduli fields generated by the string gas go like: 



where E is the energy of a constituent massless string, fj,Q is a number density and a is the scale factor of the non- 
compact dimensions. This is to be compared with the contribution of the string gas to the energy density of the 
universe, which goes like: 



We notice that the two have an inverse dependence on the energy (i.e. center of mass momentum) of the massless 
strings in question. This is a situation which would be impossible to obtain from a field theoretical framework, and as 
we will see later, quite easily allows us to evade fifth-force constraints whilst avoiding overdosing the universe without 
any fine tuning. Hence, we feel that SGC is a framework that is a priori, very well qualified to solve the moduli 
problem, the resolution of which we will turn to shortly. 
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B. Outlook, Outline and Conventions 



Since we do not wish to prejudice the applicability of this work to a particular string theory (i.e. a particular corner 
of moduli space), we leave the dimensionality of spacetime, and the number of compact dimensions general, we only 
assume that some dimensions are toroidally compactified and that F-strings (and later, D-strings) are admitted as 
stable degrees of freedom. We will, of course, pick specific cases when the need arises. The outline of the paper is as 
follows: borrowing from results in the appendix, where we derive the energy momentum tensor of a gas of strings in 
our background (following 0), we then set up the model in section 2. After carefully considering the question of frame 
switching, we study the dynamics of the coupled string gas/dilaton gravity system, and review the result that, 
given a stable dilaton, the extra dimensions are stabilized at the string scale. We then show that we obtain dilaton 
stabilization given a background with stabilized extra dimensions. Next, we demonstrate that the coupled radion- 
dilaton system also stabilizes, but in a somewhat unsatisfactory manner, in that one of the modes of the coupled 
system becomes becomes massless from the perspective of the 4-d effective theory (and hence phenomenologically 
unacceptable). We then generalize our framework to include D-string degrees of freedom, and find that in certain 
cases, possibly involving the breaking of supersymmetry, we can simultaneously stabilize the dilaton and the moduli 
describing the extra dimensions. 

Our conventions are as follows: we work in D + 1 spacetime dimensions, with the D = d + p spatial dimensions 
partitioned into d non-compact dimensions and p compact dimensions. We use the 'mostly plus' metric with signature 
( — ,+ ,+ ,..,+), and work in natural units (H = c — 1) with all dimensions expressed in terms of electron volts (eV). 
Spacetime indices will be denoted by Greek indices {(a,v...}, with the non-compact spatial indices denoted by lower 
case Latin indices from the middle of the alphabet k...}, and compact spatial indices, with lower case Latin 
indices from the beginning of the alphabet {a, 6, a..}. 



II. THE MODEL 



Consider the dilaton gravity action in the string frame (where we set the form fields to be vanishing): 




(10) 



By performing the conformal transformation: 



9^u 



(11) 



and using the relation (where tilded quantities refer to the fact that they involve g^v): 



R = e 2uJ {R + 2DV 2 lo - D(D - l)d fl ujd> M uj), 



(12) 



we can re-express the above action (after integrating the term involving V 2 



by parts) as: 



S = 




(13) 



Now, were we to set: 



(go - g) 
(D-l) ' 



(14) 



we could re-express the above action as: 




(15) 
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which is the action for Einstein gravity coupled to a scalar field. We henceforth refer to the metric g as the Einstein 
frame metric. We observe that the D dimensional Newton's constant is given by: 



8nG D = ^e 2 * . (16) 

We now introduce a gas of strings as the sources for the dilaton gravity action, which in its present form sees the 
dilaton as a scalar field minimally coupled to Einstein gravity. The energy momentum tensor for a gas of strings in 
one extra dimension was derived in and generalized to any number of extra dimensions in |nj . It was found there 
that proceeding from the Nambu-Goto action for a single string, and then hydrodynamically averaging, one would 
obtain the same energy momentum tensor for a string gas in the string frame if one were to derive this from the 
interaction term given by: 



Sint = - d D+1 x v / ~g oHoe[g]- (17) 



Where e is the energy of a single string (on which we will elaborate further on) in the given background, and /iq is 
the number density of the string gas (in the string frame) with the metric dependence factored out: 



A* = ~^=, (18) 

where g s is the spatial part of the metric. Recall that once the dilaton has stabilized, the string frame and the 
Einstein frames are equivalent. Having taken dilaton stabilization for granted, it was shown in [j| and |lfj| that a 
phenomenologically consistent stabilization of the extra dimensions at the self dual radius follows, with a resulting 
FRW late-time cosmology. 

Since it is the goal of this paper to show that, in addition to this, string gas cosmology is also capable of stabilizing 
the dilaton, we should consider the question of frame switching seriously. For this reason, and for completeness in 
general, we derive the energy momentum tensor of a gas of strings in the Einstein frame in the appendix. It turns 
out that after doing the analysis, one can obtain the energy momentum tensor of a string gas in the Einstein frame 
by transforming the above interaction lagrangian using i|ll|) : 



Sm = ~ J d D+1 x^g^e-^^ e[ge-^} , (19) 
where uj is given by (|14fl . Defining the metric of the compact directions in the string frame as 7, we see that: 



lab = e- 2u j ab , (20) 

where j a b is the metric for the compact directions in the Einstein frame. We use the result [23j that the energy of a 
closed string in the string frame, in the given background is given by: 



e = J\ Pd \2 + (n, 7 -in) + -^{w,jw) + lj[2(n,w) +4(N- 1)], (21) 
V a a 

where pd is the center of mass momentum along the d non-compact directions and ( , ) is a p-dimcnsional scalar 
product where p is the number of compact directions. n a and w a are p-dimensional vectors respectively describing 
the momentum and winding quantum numbers along the a th compact direction, and N is the string oscillator level. 
Note that in the above, we have eliminated the left moving oscillators N, using the level matching constraint @0I : 

(n,w) + N = N. (22) 

From this we find that l|19|l then becomes: 

\p d \ 2 + (n, 7 ~M + —a-(w, jw) + —7- [2(n, w) + A(N - 1)], (23) 



Sint - - / d D+1 X^ 



ffooMo 
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from which we define: 



\Pd? + (»,7 -1 ») + —=-(«>, 7«0 + -— [2(n,w) +4(7V- 1)]. (24) 
a'^ a' 



Deriving the energy momentum tensor as — -^s--£^-g\ v , we find that: 

P = ~75=, (25) 
f = ^\ Pd \ 2 /d, (26) 

p" = -^§-^}. (27) 
V5s£ b 2 a a u 

From now on, for simplicity we will drop the tildes on the metric quantites in the Einstein frame, in which we remain 
from now on. The metric tensor in the Einstein frame is taken to be diag[—l,a(t),a(t),a(t),b(t)...b(t)], d being the 
number of non-compact directions, and a' given by: 



a' = a'e 2uJ . (28) 

That is, the only effect of switching from the string frame as far as the energy momentum tensor of a string gas 
is concerned is this field dependent redefinition of a' (c.f. [T^j ). Notice that in the above, we have only derived 
the energy momentum tensor for a gas of strings with a fixed set of quantum numbers, indexed by {n, w, N, p^}, 
with a number density that will depend on these quantum numbers: \i — ^a,w,N,p 2 d - I n general our string gas will 
contain a whole host of different quantum numbers, and these will all have to be summed over to arrive at the correct 
energy momentum tensor. However it was shown in that if at some initial time we were in a situation of thermal 
equilibrium (defined through I|1H[)). then all the number densities, instead of being arbitrary are given by: 



^,N, P l=e f, **-«<"-*n r . / eP t ™t i (29) 

where the subscript 'ref refers to some reference state, which we take to have zero energy, and e is given by l|24l) - 
recalling that we have now dropped the tildes to indicate Einstein frame quantities. After careful consideration of 
what constitutes the thermal bath that couples these modes to each other, it was then shown in |lf)j | that by defining 
the following quantity (where we now generalize the result to include the contribution of the dilaton): 



ZO^A.w) :=/W ]T e -^>^W.^ (30) 

n,w, N,p 2 

one could derive the components of the energy momentum tensor as follows: 



P dp 

where P M (and in a similar way, E) is defined through the equation: 



1 3Z 1 dZ , 

Pi = T'da- = 5 ° 4 = e "' (31) 

Pa = h a §- = i§- ; &0 = e \ (32) 
P db a P d\ a 

E = 4 9 4- (33) 



^-P», (34) 
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with p M being the pressure along the fi th direction. Notice that in this setup, our interaction lagrangian (|17|l is given 
by: 



(35) 



which corresponds to the framework of Q . In and ^(j , it was argued that since the energy of a closed string in a 
given background is given by: 



' a' 



(36) 



where e is a term of order unity for the lowest lying string states, and that since for all string theories, the Hagedorn 
temperature is of order \fu' [2(J, that the summation in 13U|) projects out all but the massless string states if 
we are even an order of magnitude below the Hagedorn temperature [lpj. Since our focus is on late times, we now 
have an operational meaning of such a time- when the temperature of the universe has cooled to below an order of 
magnitude below the Hagedorn temperature. This is still considerably higher than the GUT scale, and so 'late times' 
in our framework refers only to the string scale. However this is only of passing interest to us, as we have argued that 
we shall only be interested in massless string modes from the outset. What this does indicate is that one can obtain 
these modes, demanded by phenomenological and theoretical consistency in a rather natural manner. 

Turning now to the dilaton, one finds that the associated energy momentum tensor of the dilaton field is given by: 



PD 



Pi) 



1 $ 2 



(D — 1) 8ttG d 2 ' 
4 1 $ 2 
(D — 1) 8ttGdT' 



(37) 
(38) 



where the pressure is the same for all directions (be they compact or not). Consider now the Einstein equation for 
the [i th diagonal element of the metric tensor : 



+ a M — J = 8nG D a 



Pv 



D 



D — l 



(39) 



where these equations are to be supplemented by the 00 Einstein equation: 



D D 



fl—1 v^i± 



(40) 



where is the Hubble factor associated with the fi scale factor: — a^/a^. Notice that this equation is the 
only place where the dilaton field factors, as its contribution to the energy momentum tensor drops out of l|39|) . as 
Pd = Pd- Inserting (|25() - (|27|l into (|39|l . we find for the equations of motion for the scale factors for the compact 
directions as: 



J J b=jta m 2 =0 v 



n l W 2 a b 2 a 4a, , 



e~ 4uJ e~ 2u 



b 2 a 12 



The Euler-Lagrange equation of motion for $ is given by (recall that u> = 2 ^_f^ ' 



7W) __ + w) + 2 (N 1)}] 

(41) 



2 87tGd/xo,s,^ 



— (w,jw) + —f-[{n,w) + 2{N - 1)] 

v ,jS fx 1 



(42) 



8 



The summation above is over all the massless quantum numbers which appear at the the given stage of the evolution 
of our background (recalling that (|24|l depends on the metric tensor), as previously discussed. Now we choose to begin 
at or around the self dual radius (b a ~ V«), with $ beginning at or around <&o (i-e. ui ~ as in l|14|l). We compute 
the summation over modes which are massless at, and near this point in moduli space in appendix B, where we find 
this to yield the following equations of motion: 



Ts\Pd\ 



(8p-4) 



1 



b 2 , 



-4u 



(D-r 



c=l 



-4w 



(43) 



(D-l)v^bdl 



(8p-4) 



-4w P 



-2w 



(44) 



where implicit in the above is that all the moduli are close to the values just discussed. Now we expand the above 
around: 



b a = Va'(i+r a ) ; r a «i, 



(45) 
(46) 



The equations of motion which result are: 



r o + dHT a 



3 8 ;%, ^ " 4 ) [<r. - ") - 7/7) E> - *)]' 



= 0, 



(47) 



dHu + — -4F 4(8p - 4) — 

a 3 a'va' |pd| 



a; - r c ) 



= o. 



(48) 



We now consider three separate cases of interest to us, namely when we take the dilaton to be fixed and the radii 
dynamical, when the radii are fixed, and the dilaton is dynamical, and when the radii and the dilaton are both taken 
to be dynamical. We will motivate the fact that we investiagate these cases seperately in what follows. 

A. Fixed Dilaton 

This case was already investigated in |lfj|. and so our treatment here will be brief. One might imagine that the 
dilaton is stabilized by some mechanism external to string gas cosmology, either through potentials generated by 
higher order corrections, or some other mechanism such as supersymmetry breaking. Furthermore, one assumes this 
mechanism renders the dilaton massive, with a mass that is far greater than the mass of the radion (given by the 
coefficient of T a in Q47jl). This allows us to ignore the fluctuations of the dilaton in l|47() . which results in the following 
equation of motion: 



T Q + dHT a + 



a'Va' P \pd\ 



4(8p - 4) 



r„ - 



i 



(D-l) 



c=l 



r,l = o. 



(49) 



Clearly, if we study the fluctuation of a single mode (by setting T = for all indices b not equal to the mode we're 
interested in, then the above equation describes damped oscillations about the self dual radius (r a = 0), with a mass 
squared given by the 'spring constant': 



(50) 
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However to undertake a full stability analysis, one needs to compute the eigenvalues of the Hessian matrix corre- 
sponding to this coupled radion system, and show that they are all positive. Indeed, one can write the above driving 
term as d a V , where upon defining H = d a dbV, we find that (up to an overall positive factor): 



H = X 



AD -2 -1 -1 
-1 D-2 -1 
-1 -1 D-2 

V -1 -1 -1 



where A is given by: 



- 1 \ 

-f 

-1 
D-2) 



(51) 



A = 



SttGd^ 4(8p-4) 



a 3 a'\/o 7 > £i | (D-l) ' 
Now H is a p by p dimensional matrix, which has as its eigenvalues (in units of A): 



(52) 



{D — 1,D—1, D — 1,D — p—i}, (53) 

That is, D — 1 appears with a degeneracy of p — 1, and D — p — 1 = d — 1 appears once. Clearly, these are all 
positive and hence describe a stable equilibrium. If we consider the driving term of l| l.'jll in its own right, that is by 
not restricting to small perturbations around the self dual radius, one finds that for the volume modulus, the driving 
term correpsponds to the potential indicated in fig 1. 

As found in 10], the non-compact dimensions undergo a period of FRW expansion after radion stabilization, with 
a radion mass given by the 'spring constant' in (|49() . which upon using the relation: 

G D = G d (2TtVrt) p , (54) 
and setting the scale of stringy physics to be at the Planck scale (8irG3 — 2ira'), yields a radion mass: 



™ 2 3T~7' ( 55 ) 

a 6 \pd\ 



and a 4-d energy density: 



tApd\ 

P 3~ 



(56) 
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We remind the reader that \p d \ is the center of mass momentum (energy) of these massless strings. We have neglected 
factors of order unity in deriving the above. In it was shown that for all \pd\ < lQ^eV in the present epoch, one 
can evade fifth-force constraints whilst simultaneously ensuring that the energy density of this string gas (which looks 
like hot dark matter) remains below several orders of magnitude of the critical density. 

B. Fixed Radii 

We now consider that somehow, the moduli corresponding to the extra dimensions has stabilized at the self-dual 
radius through some mechanism that doesn't couple to the dilaton (or at least contributes feebly to its dynamics), 
such as from the effects of some scalar field, or some other as yet unknown mechanism. The most obvious candidate 
mechanism, the effects higher dimensional branes in the presence of fluxes, turns out to couple to the dilaton, though 
it may be possible to arrange this coupling to not significantly affect the driving term (|44|l . The relevant point being 
that it is no less arbitrary to assume the stabilization of extra dimensions, than it is to assume the stabilization of 
the dilaton. Doing so, we find that the equation of motion for the dilation fluctuation becomes: 



8hGdHq 



(D - 1) a'^g;\p d \ 



p(8p-4) 



(57) 



which corresponds to motion in the potential (illustrated in fig 2.): 



8nG D fi 



(D - 1) a'^\p d \ 



p(8p - 4) 



(58) 



The mass generated for the dilaton by the string gas is given by the curvature of this potential at its minimum, which 
can be read off 



-4p(8p-4)- 



a 3 a' V^"\Pd\ * w {D-l) 1 
again, using (|54|l and setting the scale for stringy physics to be of the Planck scale {2na' — G3), we find that: 



(59) 



m„ 



a 3 \pd\ ' 



(60) 



whereas the 4-d energy density of the string gas that generates this potential is given (as before) by l|5t)|) . Note that 
all the features advertised in the introduction have made their appearance here. There was no need to worry about 
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exactly where the miniumum of the potentials for the various moduli occur, as there is no need to fine tune when the 
form of the action for a string gas takes the form 1|35[) . We also see that the mass of the various moduli can be made 
large enough (to avoid fifth force type interactions) without overdosing the universe (see ^3 f° r a more detailed 
treatment). In taking stock of what we have accomplished, we see that we have stabilized the dilaton around its 
expectation value (u) = — > $ = <I>o)j but we have no way of predicting precisely what this value is. Recall that the 
action (|10|l undergoes an overall rescaling under an additive redefinition of the dilaton, where the observed Newtons 
constant is then given by: 



8ttG,i 



K Q e 



(61) 



we find that whatever this $o is, the dynamics of a gas of massless strings will stabilize the dilaton around this value, 
and hence render it constant. It does not have the ability to predict what the expectation value of the dilaton will be 
(a situation that as we will see, will change when we consider the full dynamics of the system). Certainly this seems 
to indicate something peculiar about the dynamics of our moduli stabilization framework, and as we will see shortly, 
is a first indication of the presence of a flat direction in moduli space. However at present, we can view this result 
as something of a boone, in that such a mechanism is capable of stabilising the dilaton no-matter which particular 
corner of moduli space one finds themselves in. Of course one would like to be able to predict what the constants of 
nature are from a given moduli stabilization mechanism, and as such this is somewhat of an unsatisfactory situation. 
It is only of passing interest to us however, as the main focus of this report will be the case of when we consider the 
coupled dynamics of the radion-dilaton system, to which we turn our attention to now. 

III. COUPLED DILATON- RADION DYNAMICS 

Recall the equations of motion for the fluctuations of the radion and the dilaton fields l|47l) (|48l) : 



f a +dHT a + 8lTG ^ r 4(8p-4)[(r g -o;)- ' 



p 



(r c - w)] 



= o, 



(62) 



LO - 



dHCo 



8ttG 



4(8p - 4) 



(D 



hi 



v 

E 

i 



(w - r c ) 



(63) 



We have already see how the radion fields and the dilaton are separately stable. It is a straightforward generalization 
to study the case of when both sets of moduli are dynamical. To find out if the coupled system is stable, we have to 
find out if we are at a point in moduli space which corresponds to a global minimum. It seems straight forward to 
generalize the treatment above, where it seems that all we have to do is to compute the eigenvalues of the Hessian 
matrix. However there is a subtle issue here- the driving term for the dilaton is given by: 



dZ_ 

doj ' 

whereas the driving term for the radii is given by (c.f. lj3U|) - (|33)l ): 



(64) 



D - 1 



a o \ 
^— ' v da v 36 J 



Z. 



(65) 



If one were to write the driving term for the radii as 4p for some V, then in general this will not be the same V that 
one derives the driving term for u> because of their different origins (wheras one comes from linear combinations of 
the Einstein equations, the other comes from the Euler Lagrange equations). However, were one to make the trivial 
rescaling: 



(66) 
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where d is the number of spatial dimensions, one can write the above driving terms as the partial derivative of a 
'master potential'. The neccesity in doing this comes from the fact that were we to naively take the second derivatives 
of the driving term above, one finds that d u dr a V ^ dr a d LO V. This rescaling allows us to write the above driving 
terms as the derivative of the appropriate potential. From this, we find that, up to an overall positive factor: 



d a a u v = d u d a v = -j§~Yy ( 67 ) 

BlV = jJL-y (68) 

d a d b V = -j^Yy (69) 

9 2y ( £> ~ 2 ) (70) 

which corresponds to the p + 1 by p + 1 dimensional Hessian matrix (again up to the positive factor l|52|0: 



H F = A 



( D-2 

-1 
-1 



-1 
D-2 
-1 

-1 



-1 
-1 
D-2 

-1 



\-Vd~l -Vd~T -VcT^l 



-1 - Vd~l \ 
-1 - Vd~T 
-1 -Vd^T 



D-2 -y/d~T 
-Vd~T p j 



(71) 



This matrix has eigenvalues X(D — 1), which appears with degeneracy of p, and 0, which appears once (in fact the 
corresponding eigenvalue for such a p+l-dimensional matrix is X(—p — d + D) = 0). The appearance of a zero mode 
to this system implies that one degree of freedom in this coupled system is not stabilized (albeit marginally). In fact, 
one can directly obtain the normal modes of the system from i|t)2ll 1631) : 



K + dHA a + - -£^4(8p-4)A a = ; A a :=T a -u, (72) 

a d aWoi' \pd\ 

p 

k + dHk = ; k = ^2r a + (d- 1)uj, (73) 

a=l 

from which we conclude that the individual radial moduli (r a ) tend towards the dilaton (i.e. T a — lu —t 0) quite 
robustly, which, from the definition of k, results in k tending towards (D — 1)lu, which now no longer has a mass. 
This is not to say that this modulus is free to run away, as according to (|73|l . one can solve for n to find that: 



«(t) = c/J^ (74) 

which implies that in a universe with expanding non-compact dimensions (which is easily obtained in our framework 
this modulus terminally approaches some finite value. However, this is far from satisfying to say the least, 
as this modulus is no longer robustly held at any particular value and fluctuations and interactions are likely to shift 
the value of this modulus significantly through time. 

One can quite easily understand the physical origins of this zero mode. Recalling the energy of a closed string l|24(l 
in the given background: 



e = \ \p d \ 2 + (n,7- x n) + — ^-{w,jw) + — — [2(n,w) +4(N- 1)], 
V a a 

we see that if instead of stablizing at 0, were u> to stabilize at some other value, say wo, then modes which were 
massless for uj = at the self dual radius are now massless at u> = ujq at the radius defined by b a = \fa! e"° . Hence 
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were we to begin at this point in moduli space, and compute the sum over massless modes in our equations of motion 
(as in how we arrived at (|62[H|()3[l ). one would end up formally identical equations for the fluctations about this point 
in moduli space, and obtain the same results for our stabilizty analysis. This would certainly indicate that we have a 
flat direction in the dynamics of our coupled system. 

Before we attempt to resolve this unsatisfactory state of affairs, we wish to point out that were the dilaton to 
receive a potential due to effects external to SGC (i.e. through higher order corrections), then one retrieves a stable 
equilibrium for our system for any value of the stabilized dilaton to = loq, where the extra dimensions stabilize at 
b a = \fa' ^e u ° . That is, were some external potential to stabilize the dilaton, the Hessian matrix at the stabilized point 
in moduli space would then read: 



/ 


D - 2 


-1 


-1 


-1 


-Vd-1 \ 




-1 


D -2 


-1 


-1 


-Vd-1 




-1 


-1 


D-2 .. 


-1 


-Vd-1 




-1 


-1 


-1 


. D-2 


-Vd-1 


V 


-V<2-1 


-Vd-1 


-V<2- 1 •• 


. -V<2-1 


P + cV"{lu)J 



where the only difference is in the last diagonal entry. In the above, c is some positive factor (the overall factor for 
the string gas contributions to the Hessian will not be the same as that for the dilaton stabilization mechanism H52|) ~). 
This matrix generically has all of its eigenvalues as positive for any positive V" . For instance, for 6 extra dimensions 
in D = 10, one finds that the Hessian matrix has the eigenvalues (in units of A): 



{9, 9, 9, 9, 9, -(9 + cV" - a/81 + 6cV" + (cV") 2 ), -(9 + cV" + y/81 + 6cV" + (cV") 2 )}, (76) 

which are all positive. Hence we conclude that one can consistently stabilize extra dimensions in string gas cosmology 
independent of any particular (external) dilaton stabilization mechanism. This ameliorates the difficulties encountered 
in [U, which studied precisely this aspect of the problem, but neglected the effects of massless string states. 

Returning however to the issue at hand, which is to explore ways in which some mechanism within SGC might 
lift the zero mode of the coupled dilaton-radion system, we offer a quick recap of the possibilities raised previously. 
Firstly, one could introduce higher dimensional branes and turn on fluxes, in the hope that this can somehow enhance 
the stabilization of the extra dimensional moduli but not the dilaton (or vice versa), in which case we would end up 
with the two separate situations considered in the previous section. Alternatively, one could introduce branes and 
fluxes and not aim for preferential enhancement of the stabilization mechanisms of one of the two sets of moduli, but 
instead see if doing so can somehow otherwise lift the zero mode (as we have seen, obtaining even a small additional 
positive component to the last diagonal entry in (|75(l suffices). However there are compelling reasons for us to focus 
solely on one particular possibility- namely the introduction of a gas of D-strings. 



IV. D-STRINGS 



In the introduction, we made a rather detailed case as to why massless string modes should be the only concern of 
string gas cosmology if it aspires to be relevant to the present epoch. Consider then what it would mean to introduce 
branes of various dimensions into our framework. Since we do not yet know how to quantize any extended object 
other than strings, we will have to content ourselves with a purely classical treatment of such objects, for which the 
lowest energy configurations have energies: 



E = T p x Vol (77) 

where Vol is the worldvolume of such an object, and T p is the tension of such a p-dimcnsional object (typically the 
p'th power of the string tension). We see that introducing a gas of such objects would run counter to the spirit of 
what we have done so far, in that we would then be introducing objects which have a significant energy density, calling 
into question any applicability to late time cosmology, and possibly taking us away from the regime of the low energy 
effective theory (|15f) in which we are working in. Hence if we are to be consistent with our own demands, we should 
ignore these objects altogether. However, we find that there is one class of D-branes which evades the objections just 
made, namely D-strings. Not only are D-strings quantizable (in certain backgrounds, with caveats), they also have 
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the same massless excitations as the F-string |23j , and offer a window of opportunity to lift the zero mode uncovered 
in the previous section as they couple to the dilaton in a different manner than F-strings. This follows from the action 
for D-string in the absence of fluxes (in the string frame): 

S = -^J d 2 ae-*V^h (78) 

where Ti/2ir is related to the D-string tension (:= t\) in a background with the dilaton assuming the value $o through 
the relation: 

ri = (79) 

Z7T 

We will elaborate on the tension of the D-string further on. In what follows, we will be using the fact that as mentioned, 
one can indeed quantize the D-string in certain backgrounds. As a solitonic object in string theory, whether one can 
do this or not is a highly non-trivial question. However it should be clear that for a constant dilatonic background, 
(I78|l describes the same action as that for the F-string, with a corresponding redefinition of a' . One should certainly 
be able to quantize this action and recover the spectrum of states (which follows from the constraint algebra alone), 
if one can do this for the F-string. It was shown in [25l ] that one can indeed covariantly quantize the superstring, 
provided the form fields satisfied certain conditions. In [2(j the canonical equivalence of the constraint algebra of 
the D-string and the type IIB string is established, which implies that these two objects have the same spectrum. 
Since we have endeavoured to keep the background theory in which our framework is established arbitrary, we do not 
consider this issue further, taking as an assumption that one can indeed quantize D-strings in the corner of moduli 
space we are considering. We argue in the appendix that for situations in which the dilaton field changes much, much 
slower than the string scale (defined by <E> <C l/2na'- which we show later holds to an excellent approximation in 
our model), the only change in the D-string gas energy momentum tensor after switching to the Einstein frame i|ll|l . 
compared to that for the F-string gas H25|) - (|28|) . is the replacement a' — > e*/Ti: 

P = (80) 
p* = -^IPdWd, (81) 



PD ,n\ 



rgl^bi (e*/T x ) 



2 



(82) 



(83) 



Further issues concerning the consistent introduction of a D-string gas are considered in section V. Recalling the 
relation i|14l) . one can rewrite $ as: 



* = $o - -— ^> ( 84 ) 



and by defining the following quantity: 



a-.^e^/Tx, (85) 
we can rewrite the above expression for the energy of a single string as: 



e (D-5)u e (D-5)w/2 
— (W,7W)H — 

a z a 



( = \l\Pd\ 2 + (n,7 _1 n) + — — Ktw) + [2(n,ifl)+4(JV-l)]. (86) 
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where the numerologically peculiar factor D — 5 will turn out to have little consequence. A gas of D-strings (with 
a fixed set of quantum numbers) will yield an energy momentum tensor (in the Einstein frame) that can be derived 
from the action: 



f n,, / e (D-5)u e (D-6)u>/2 

S = - d D+1 x^p D ^^ N \\ Pd \ 2 + (n, 7 - 1 n)+ ~ 2 (w,jw) + 2(n, w) + 4(N - 1) , (87) 

7 V or a 

where the subscript D on p, indicates that this number density is specific to the D-string gas, which for massless states 
we take as equal to the number density for the F-string gas. Just as before, if we realise that one has contributions from 
each massless state, then one has to perform a sum over massless modes to arrive at the correct energy momentum 
tensor (see discussion immediately following J25J)). The equations of motion that result are as follows: 



8irGDfJ-D,n,wb a 



b 2 a a 2 + (U-1) LW7W; 



(D-5) 



,(D-5)w/2 



-{(n,w) + 2(N-l)}} 



(D - 5) 87rGDMD,n,t« 



^ 2(D 

m 2 =0 y 



1) 



(D-5)u>/2 



[(n,w)+2(N-l)} 



(89) 



We now make the observation that at a generic point in moduli space, quantum numbers which render the F-string 
massless are not the same as those which render the D-string massless (a detailed treatment of what follows is given 
in appendix B). For instance, at the self-dual radius (7 & = a'6 a i,), with the dilaton stabilized around <& (i- e - w = 0), 
the mass of the F-string is given by (c.f. J2U0 : 



a'rn 2 = (n + w,n + w) + 4(N - 1), (90) 

from which we catalogue the classes of states which are massless (recalling that the above formula must be supple- 
mented with the level matching constraint N + (n, w) > 0): 



N 


(n, n) 


(w, w) 


(n, w) 


1 











1 


1 


1 


-1 





1 


1 


1 





1 


3 








3 


1 








2 


2 








4 














4 






Whereas the mass of the D-string at this point in moduli space is given by: 

a'rn 2 = ( n ,n) + 2 (w, w) + 2^) [(n, w) + 2{N ~ 1)], (91) 

where a is defined in 1)85(1. We see from this that unless a = a' , the quantum numbers tabulated above are not going 
to correspond to massless D-strings. Since the stabilization mechanism uncovered previously |10| relied so heavily on 
the features of modes with these quantum numbers, right down to their degeneracies and behaviour away from the 
self-dual point (see appendix B), we take as a non-trivial assumption that we are in a region in moduli space where 
indeed a — a'. Doing so allows us to sum over the same set of quantum numbers in (|88|l and (|89|l . and as we are 
about to see shortly, this will have some striking consequences. We will return to address the non-trivial point of 
whether or not one can consistently set a to be equal to a 1 further on. For now, we see that the effect of introducing 
a gas of D-strings at this particular point in moduli space is to contribute to the equations of motion (|88|l(|89|l as: 



b a 



iTrG D n D b a 
V<h\Pd\ 



(8j>-4) 



1 h 2 
b 2 , a 12 



,(D-5)u>/2 



(92) 
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(J - 5) 8nG Df i D 
'2(D-1) VSHPdl 



(8p - 4) 



rp (B-5)w P 



,(D-5)«;/2 



(93) 



Of course, implicit in the above is that all the moduli are close to the values just discussed (all radii are close to the 
self-dual radius and $ ~ consisitent with a — a'). Just as before, we expand the above around: 



b a = vv(i + r ) ; r Q «i, 

to ; 



The resulting equations of motion are: 



(94) 
(95) 



c=l 



= 0, 



(96) 



By performing the same rescaling as in (|66|l . one can write the driving terms of this coupled system of equations as 
the derivatives of an appropriate potential, from which we arrive at the Hessian matrix (up to the same positive factor 
(|52|l that appears in front of J7TJ): 



/ D-2 
-1 
-1 



H D = \ 



-1 
D-2 
-1 



-1 
-1 
D-2 



K^Vd^i ^Vd^i ^Vd^J 



-1 
-1 
-1 

D-2 



(98) 



Hi 



/ 



However, we have stipulated that both F-strings and D-strings will be present, so in addition to the driving terms 
generated by the D-string gas in the above, we will also have the driving terms generated by the F-string gas l|47|l . 
(|48(l , which will result in the Hessian matrix for the combined system to be Hd+ f = Hf + Hp , where Hf is given by 
H7ip. It is here that we find the result we have been looking for: the eigenvalues of Ho + p are all positive, inspite of 
both Hp and Hp individually possesing a zero mode. Hence all moduli fields have a stable equilibrium at the point 
in moduli space we have identified. For the specific case of D = 9, p = 6, we find that the eigenvalues of the Hessian 
matrices we consider are: 



H F : {8,8,8,8,8,8,0} (99) 
H D : {8,8,8,8,8,8,0} (100) 
H F+D : {16,16,16,16,16,12,4}, (101) 

where all eigenvalues are in units of the overall positive factor (|52|l . We wish to point out that the fact that Hp and 
Hp, have identical spectra is only a numerological coincindence for D — 9,p = 6. One can investigate the stability of 
the coupled D-string and F-string system in any number of dimensions (albeit numerically), and find that we always 
lift the zero mode we uncovered previously. For instance, for D = 25, p = 22, we find: 



Hp : {24,24,24, ...,24,24,0} 
Hp : {552, 24, 24, ...,24, 24,0} 
H F+ p : {12(24 + \/565),48,48, ...,48,48, 12(24- \/565)} 



(102) 
(103) 
(104) 
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Where the last eigenvalue in Hp + £, evaluates at 2.76A. What this implies is that all of the moduli we consider, are 
stabilized at the self-dual radius (b a = Vol), with the dilaton stabilized at the point determined by the condition 
a = a'. Before we consider the phenomenology of this conclusion, we should return to address the issue of whether 
or not one can consistently set a' = a. 

We know in general that the D-string tension is related to the F-string tension (= l/2ira') through the following 
relation, which is exact and is a consequence of supersymmetry |23| : 

Tl = w" (105) 

Therefore, one can only set the two tensions to be equal to each other at the self S-dual point (e*° = 1), from which 
we conclude that at weak coupling (which is the only regime in which the low energy effective theory is valid), 
it is impossible for the two tensions to be equal to each other. Hence the stabilization mechanism we just proposed 
would seem to be untenable. In fact, for a generic point in moduli space corresponding to weak coupling, we have the 
relation: 



%■ > 1. (106) 
a 

Thus one will never have a situation in which the same quantum numbers describe massless excitations of both the 
wound D-string and the wound F-string. It is interesting however to determine what the effects of the massless 
states of the D-string are, at a generic point in moduli space consistent with the weak copuling limit. Clearly states 
with the quantum numbers N — l,n,w = correspond to massless D-strings as they do to massless F-strings, and 
correspond to unwound gravitons, which do not directly contribute to the dynamics of the radion fields or the dialton 
(c.f. ^88 p .(|89 p for the D-string and (|41|l . i|42[) for the F-string). If were we to start out at the self-dual radius at 
weak coupling, then in general the only D-string states which are massless correspond to unwound graviton states 
N = 1, ft, w = 0. However, for the specific case when a' jo. — q; q e Z ^> 1, it can be shown that the only massless 
states are those with quantum numbers: 



(n,n)=4q; N = 0, w = 0, (107) 

the specific realisations (and degeneracies) of which, depend on the number of extra dimensions we are considering. 
We see that quite generally, these states consist of only momentum modes and from (|88|l and 1)89)). we see that they 
have the net effect of causing the radion fields and the dilaton to run away. The run away behaviour for the dilaton 
field is such that it drives the dilaton towards stronger coupling. In particular, if we start out at weak coupling, the 
dynamics of the system tend to drive it towards the self dual coupling point- where one would naively expect to find 
moduli stabilization from the results just uncovered. Quite obviously however, one expects qualitatively new physics 
to appear at the self S-dual point in moduli space, and so our framework is expected to break down at this point as 
it is highly unlikely that these new effects will conspire in such a way that it leaves our framework alone. 

Since the issue of the D-string tension i|l(J5[) is not up for debate, it might seem as if we have failed in our efforts 
to find a convincing moduli stabilization mechanism withing the framework of SGC. However, we have to keep in 
perspective that not only are we attempting to solve the moduli problem, we are attempting to solve it in a way that 
is relevant to the present epoch. Specifically, we are proposing a model which to our knowledge is the only attempt 
at stringy cosmology which makes contact with, and is operational in the present stage in our universe's evolution. 
We made a case for why one should only consider massless string degrees of freedom in the introduction with this in 
mind, and with this in mind presently, one should certainly expect us to be in a situation where we are in a regime 
with broken supersymmetry. It is almost immediately built into our framework that we are in a situation of broken 
supersymmetry, as introducing a gas of massless D-strings, whose constituents are moving around at the speed of 
light with respect to each other, are not going to be BPS (many bodied) states. There are exactly two consequences 
of this observation on what we have done so far. Firstly, since we have restricted ourselves to the bosonic sector of 
whatever string theory one might wish to apply this model to, the only difference for our framework will be in the 
critical dimension in which we wrote down the low energy effective action Since we have deliberately left this 
dimensionality arbitrary, we do not have to change any aspect of our setup. However, one can make the argument that 
one can still consider spacetimes with the usual dimensionality (i.e. either 10 or 11) in spite of broken supersymmetry, 
for if we can indeed show stabilization at the string scale that is consistent with the low energy framework, any external 
physical probe of the scale of these extra dimensions will neccesarily have to be of a very high energy. At such energies, 
one certainly expects supersymmetry to be restored. This allows us to conclude that at most, the critical number of 
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dimensions dictated by superstring theory are observable, inspite of the present epoch having broken supersymmetry. 
However this last point is in no way essential to our argument, as we have shown moduli stabilization in any number 
of spatial dimensions (D=10, 11, 26...), provided one can find a point in moduli space where the D-string and the 
F-string tension are equal to each other. 

The second consequence of broken supersymmetry is by far the most important to us, for in this case, the relation 
(|105|l is longer true. In its place, is the relation (derived in D=25, by comparing F-string exchange amplitudes between 
two D-strings) [^ p i f : 

ti = 1fi ^^ ttV) 5 , (108) 

where kq^d is the normalization of the D-dimensional low energy effective action (J2J, which is determined by studying 
dualities between various string theories, and is proportional to a' to the appropriate power. The observed D- 
dimensional Newton's constant is given by: 

8itG d = Ko. D e 2 * . (109) 

If one were to have compactified (or achieved radial moduli stabilization of) the p extra dimensions, kq satisfies the 
relation: 

<D = nl d a'^VolT p , (110) 

where T p is the volume of the p compactified dimensions in units of a' p / 2 (T p — {2tt) p for toroidal compactifications) , 
and the observed d-dimensional Newtons constant (d being the number of non-compact dimensions, which we will 
soon take to be 3) is given by: 

STrGd = 4 d e 2 * . (Ill) 

We immediately see from substituting (|110|l into i|108|) . that the D-string tension depends on the compactification 
mechanism we have through the new relation: 

Vi(47rV)5 

1 16K ,da'f/ 4 (^T p ) 1 /2 e *o ' ^ ' 

which indicates the possibility of tuning the tension of the D-string, which depends on the string coupling, the 
compactification scheme and the particular string theory we are working in. Returning to the original expression 
(|108fl . let us now assume that we start out at a value for $o, such that the D-string tension is equal to the F-string 
tension: 

n = J- -> e*o = f^(4,V) 5 (113) 
lira' lbfto^ 

The overall coupling kq.d that appears in the low energy effective action J5J) is given for type I, IIA, IIB and Heterotic 
string theories (in D = 9) as [23| : 

K§ )9 = i(27T)V 4 . (114) 

We thus conclude from the general relationship i|110|) , that when we are in a situation where all extra dimensions are 
toroidally compactified at the self dual radius, that the relation (|108fl allows us to abstract from i|113|) the value of 
the string coupling where the D-string and the F-string tension are equal to each other: 

e*° = 1. (115) 
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In this case, one would have the situation where oi = a consistent with weak coupling, and massless excitations of 
the D-string would have the same quantum numbers at the self-dual radius as those of the F-string. We conclude 
that in this case, it is certainly possible for a gas of massless strings to stabilize all the moduli we have considered. In 
fact, as argued earlier, if we were to begin away from this point in moduli space, the dynamics of the coupled dilaton 
gravity-string gas system would drive the moduli fields to this point, where stabilization would result. We also wish 
to take note of the fact that the usual instability of D-strings in bosonic string theory is no longer at issue when the 
tensions of the two objects are equal, as the distinction between which one is fundamental becomes blurred, as there 
are no longer any obvious decay channels for the D-string. Note that not only does this mechanism make a prediction 
for the string coupling constant: 



9s = ^, (116) 

it also predicts the observed 3-dimensional Newtons constant (through the appropriate dimensional reduction of l|114|) l 
in terms of the string scale: 



9-77- <V 

8-G3 = (117) 

which is certainly in the right ball park if we set the scale of stringy physics to be at or around the Planck scale. The 
phenomenology of this stabilization is exactly the same as the special cases we considered in the previous sections. 
All of the moduli fields, when expressed in terms of normal modes will have masses given by the curvature of the 
stabilizing potential at the minimum, given by i|52|) , which upon using the relation (|54|) gives us masses for the moduli 
as: 



2 ^0 / 1lo i 
m ~ — (118) 

a 6 \pd\ 



Similarly, the contribution of the string gas to the energy momentum tensor will go like: 



P 3— , (119) 

from which it was shown in || and ^(Jj that for \pd\ < 10~ 3 eF, it is possible to evade fifth force constraints whilst 
simultaneously ensuring that the energy density of the string gas is substantially less than the critical density of the 
universe. Notice that the only parameter we have to tune this model is the center of mass energy of these strings 
(IpdDj and even this only has to satisfy a rather mild bound than to be tuned in any sense. Everything else in this 
framework determined by the requirement that we are in a regime where low energy effective theory is operative, and 
hence only massless modes will be contributing to the dynamics of spacetime. It is through the manifold remarkable 
properties of massless modes that we obtained our desired result. We wish to make the remark that string modes 
which have quantum numbers associated with massive modes not only have an inconsistent phenomenology, but 
are also not likely to be around in any significant number in the present epoch, in addition to not leading to any 
stabilization mechanism. This is because any potential generated by massive string states refuses to admit terms 
which compete with one another such that a minimum results. It even appears that a no-go theorem to this effect, 
for moduli stabilization using massive modes is afoot |27| . However, pressing on with our goal of stabilizing all the 
moduli of our given compactifications, we now turn our attention to stabilizing the shape moduli for the toroidal 
geometries we've been examining. Along the way we'll tie up a few loose ends, after which we'll conclude this report. 



V. SHAPE STABILIZATION, CONSISTENCY CHECKS AND SUMMARY 



Since the issue of stabilizing the shape moduli was already considered in |13( , we do not consider the issue in great 
detail here as the mechanism through which it stabilizes is identical in our setup. In general, when we refer to the 
shape moduli of the toroidal compactification we are studying here, we mean the off-diagonal components of the 
metric for the internal manifold. Allowing the shape moduli to be dynamical means that our toroidal metric has 
components such that: 



lab ^ , Va^b. 



(120) 
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For illustrative purposes, we consider the simple example of a compactification on a two torus. How the results that 
follow generalize to higher dimensional tori will be indicated as we progress. In this case, the metric of the torus has 
the form: 




(121) 



If we consider turning on the shape moduli such that they are close to vanishing, we find that to first order these 
moduli fields only make a difference to the off-diagonal elements of the Einstein Tensor (a result that is also valid 
for higher dimensional tori). Hence all of the results just derived for the other moduli fields are still true- the shape 
moduli decouple completely from all other moduli to first order. In the case we are considering, we only have one 
extra Einstein equation to consider, namely the off-diagonal equation for the internal dimensions. This equation will 
read: 



c+(dH-Hi-H 2 )c+ [c(J2 P M - P) - 4Tia] = 0, (122) 

M 

where we picked the indices 1,2 for the compact directions. The term (X^^Pm — P) vam shes around the moduli 
stabilization point, and deviations from this point result in terms which are linear in the infinitesimal fluctuations of 
the radial and dilaton moduli fields. Since c itself is infinitesimal, this term is negligable. We derive T\ 2 for a gas of 
strings with a fixed set of quantum numbers as: 



Tn = -c- 



/'-() 



-4w 



(123) 



Again, since c is infinitesimal, we conclude that the above term contributes to <|122|) as: 



c+idH-Hx- H 2 )c + c 8 -^ = 0, (124) 

where N is some positive number which corresponds to the degeneracy of such a term once we sum over all massless 
modes. Again we see that this modulus is stabilized. Using 1)54(1 . and setting the scale of stringy physics to the Planck 
scale, we compute the mass of this modulus to be: 



2 PO /, rl 

m ~^' (125) 

which is similar to the masses found for all the other moduli fields (I118(l . The case for extra shape moduli, such 
as what would appear for higher dimensional tori, should be no different, but nevertheless should be examined. We 
postpone this study to a future paper. 

We conclude by addressing a few points concerning the consistency of this framework alluded to in the introduction, 
and offering our closing thoughts. The first thing to check is that we have remained within the regime of low energy 
effective theory throughout the evolution of our system on the way to moduli stabilization. This is easily done by 
taking the trace of Einstein's equations, from which we find that: 



n = ^t;. (126) 



Refering to 1(39(1 and 1)43(1. we see that this implies that when the system is fluctuating around the stabilization point 
in moduli space (again using ijSU) . and setting the string scale to be near the Planck scale): 



(127) 
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where the fluctuations of the moduli fields will appear in the parentheses. Clearly, since the prefactor is of the order 
of electron volts (so that we evade fifth force constraints- see (I118|l 9] 10] ), one quite easily satisfies the condition |J3}: 

i?<^r. (128) 
a' 

We make the comment here that were we to consider massive string modes as the source terms for the dynamics, one 
typically could not satisfy i|128|) . 

Another aspect of our model that should be checked, is that the energy momentum tensor is covariantly conserved: 

D 

V»T<? = p + ]T #„(p + p„) = 0. (129) 

It is easy to check that this is (quite remarkably) satisfied as an identity by the coupled string gas-dilaton system 
(|25 [l -l]27 [l and (|37fl - i|38|) . In taking the time derivative of p, which is a function of the metric and the dilaton, one 
finds that the time derivatives of the metric factors are cancelled by contributions from the terms proportional to 
the Hubble factors in the above, and the terms involving time derivatives of the dilaton are cancelled by the same, 
subject to the equations of motion for the dilaton field. 

We also wish to note that since we have chosen a simple toroidal compactification for our compact manifold, one 
does not have any external sources for the {B^ v ) flux fields, as would be the case if we were to have used an orbifold 
compactification (the orbifold fixed planes would then source the flux fields). Hence it is consistent for us to set these 
fields to be vanishing. A further consequence of this compactification scheme, is that a gas of D-strings had better 
contain an equal number of RR charges and anti-charges if Gauss's law is to be satisfied. Since we are considering 
low energy phenomena, any hydrodynamical fluid of such a collection of D-strings will look RR charge neutral, and 
will not source RR fields, which also allows us to set these fields to be zero. 

As this concludes the body of our report, we feel that a summary would be in order at this point to realign our 
perspective before we offer our closing thoughts. We have demonstrated in this paper, that it is possible to stabilize 
all of the moduli in the context of a general toroidal compactification in the context of bosonic string theory. The 
workhorse of this stabilization mechanism is a gas of massless F-string and D-string modes. In the introduction, we 
provided several a priori arguments as to why one should focus exclusively on these modes, after which we derived 
the energy momentum tensor of a general string gas. We then showed how one quite naturally ends up with only the 
massless string states if we are at an energy that is even one order of magnitude below the Hagedorn temperature. Upon 
studying the dynamics of the various moduli fields coupled to this string gas, we found that we obtain stabilization of 
the radial moduli, given the stabilization of the dilaton. We then found that we obtained dilaton stabilization, given 
the stabilization of the radial moduli. However, upon examining the dynamics of the coupled radion-dilaton-string 
gas system, we found that a zero mode appears for one of the moduli. The physical origin of this zero mode was 
uncovered, and we then looked for ways of lifting this zero mode. We used the arguments presented in the introduction 
to focus on one possibility- the introduction of a gas of D-strings. We found that at the point in moduli space where 
the D-string tension was equal to the F-string tension, stabilization of all moduli resulted. We then examined the 
question of how one would obtain a situation where the tensions of the two classes of strings were equal, and found 
that in a situation of broken supersymmetry (as would neccesarily be the case in the present epoch, especially if we 
had a gas of massless D-strings propagating through our space), that this was certainly possible consistent with weak 
coupling. We then considered the consistency of our framework with the low energy effective theory that provides 
us with our gravitational theory, and found that our model passes all of the tests we put to it (phenomenological or 
theoretical) quite easily. Hence, albeit in the context of a specific model, we feel that we have presented a model of 
late time cosmology that offers a solution to the moduli problem. We conlude with a discussion of the strengths and 
weaknesses of this approach, and open questions that this new mechanism begs. 

VI. CONCLUSIONS AND OUTLOOK 

In the course of this study, some rather remarkable properties of massless string modes were uncovered. Not only 
did these properties prove crucial in stabilizing all of the moduli in our framework, but had, in addition, many other 
desirable phenomenological features associated with them. In general, we expect plenty of such strings to be present 
in the present epoch. As massless objects, these are likely to be copiously produced even if a preceding epoch of 
inflation took place. A very desirable feature of our framework is the distinct lack of fine tuning involved. Unlike 
most other attempts at moduli stabilization, which require fluxes to be turned on in a certain manner, with specific 
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brane configurations to be present, introducing a string gas allows no such freedom. From the form of the coupling 
of the string gas to gravity ©, we see that we do not need to tune the minima of our 'moduli potentials' to avoid 
an unacceptably large De-Sitter phase of evolution. Another advantage is that evading fifth-force constraints without 
overdosing the universe is easily accomplished by a string gas, again without any fine tuning (see (|118|l . (|119|l . 
and [10(). In [1 1| . we will investigate further the observational consequences of having an all pervasive gas of strings 
stabilizing extra dimensions at the present epoch, which we feel might be within the reach of experiments in the near 
future. However all of these results are just a desirable aside to the main result, that we can consistently stabilize all 
of the moduli fields in our compactification within our framework. 

It is prudent at the very least, to address the shortcomings of such an approach, not only to inform future avenues of 
investigation, but to also contrast to other approaches to moduli stabilization (see [| and references therein) . Clearly, 
one cannot hope to stabilize moduli fields associated with the compactifications of string theory solely using string 
gases if the the compact manifold that results does not admit topological 1-cycles (see however 0). That is, since 
1-dimensional extended objects need 1-cycles to admit winding quantum numbers, manifolds that do not admit these 
(such as Calabi-Yau spaces) will not readily admit to our mechanism. In addition, we have not yet presented any set- 
up that generates an inflationary epoch for our non-compact spacetime within our moduli stabilization mechanism 
(however there are promising indications one such a mechanism exists |12I . see also |17|V We presume that this 
situation will be markedly different if our framework was generalized to include brane gases. For in this case, not only 
can we hope to stabilize compact sub-manifolds that do not admit 1-cycles (but do admit higher dimensional cycles) 
through the use of higher dimensional branes, but one might also expect inflation to fall out of such a set-up. One 
can envisage that the various inter-brane potentials that arise might generate an inflationary epoch on the way to 
moduili stabilization. However in this case, the framework we propose starts to look a lot more like the more standard 
approaches to moduli stabilization Q , and it is highly likely that these approaches will overlap significantly. It is our 
belief that the mechanism that will eventually succeed as the main candidate for stabilizing all of the moduli fields 
of string theory, will certainly be some combination of the flux-compactification and the string/brane gas approach. 
However, insofar as providing a self contained, consistent and natural mechanism for moduli stabilization in the 
context we have chosen to study, we feel that the string gas approach has several comparative advantages over the 
flux compactification approach. It is our hope that, at the very least, the results uncovered in this report will help 
motivate the further study of string gas cosmology, not only as a mechanism to stabilize moduli, but also as a viable, 
alternative stringy cosmology. 
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APPENDIX A: THE ENERGY MOMENTUM TENSOR 

We offer in this appendix, a derivation of the string gas energy momentum tensor based on micro-physical consid- 
erations. That is, we derive the space-time energy momentum tensor of a single string in a given background, and 
proceed to hydrodynamically average to obtain the result for a gas of such strings. Before we do this, we should 
comment on the fact that we are going to be using aspects of the string spectrum in a background that is varying 
in time. Certainly, string theory in a time-dependent background is an active and current area of research, in which 
many unresolved issues remain. However, we are going to be doing string theory in which the background varies on 
a time scale that is many, many orders of magnitude slower than the string scale, and as such one would expect to 
inherit many features of the string spectrum in this very weakly time dependent setting. This intuition is justified 
rather more rigorously in where we discuss quantization of the string in such a time dependent background and 
show that (in what is essentially an application of the adiabatic theorem of quantum mechanics), we have the same 
spectrum of states in a weakly time dependent background, with the energies now very slowly varying (compared to 
the string scale) functions of time. We arrive at this conclusion from considering how the constraint algebra of string 
theory changes in the context of a weakly time-dependent metric. 
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Returning to the problem at hand, we begin with the Nambu-Goto action: 



S = -^JM~K (Al) 
where the induced worldsheet metric h a b is given by: 

h ab = g ia ,{X)d a X' l d b X v , (A2) 
and the spacetime metric (in general a function of the worldsheet fields X 1 - 1 ), is assumed to have the form: 

.900 = ,9oo(A°) (A3) 
9ij = <%a 2 (X°) (A4) 

gab = lab(X°). (A5) 

By varying (|A1|) with respect to the background metric, one obtains the spacetime energy momentum tensor as: 

9 

TP" = Z_ r ■ (A6) 
An arbitrary variation of the background metric induces a variation of the worldsheet metric as: 

S Xf, g ltv = 6*558 D+1 (X T -y T ) - 5 xr3 h ab = d a X x d b X' 3 5 D+1 {X T - r /), (A7) 

where the unmatched indices A and (3 imply that we perturb only these components of the metric tensor, and 
S D+1 (X T — y T ) is a delta function in D + 1 spacetime dimensions. The result of this variation yields an expression 
for the energy-momentum tensor for a single string as: 

T"" = { d 2 a^hh ah d a X^dbX v 8 D+x {X T - y T ). (A8) 

^P^gl-Ka' J 

Before we can proceed, we must pick a gauge to work in. We choose to work in conformal gauge, defined by: 

hab = A ^ JJ , (A9) 

where we keep this up to an arbitrary positive factor, which we will choose further on. From (|A2|) . we see that this 
gauge choice implies the conditions: 

g^X^X" + 9llv X'»X ,v = (A10) 

g^X^X* = 0. (All) 

It can be shown that even though we are in a (weakly) time-dependent background (in particular one that is not flat), 
one can still make this gauge choice simultaneous with the condition: 

X'° = 0, (A12) 

where the prime denotes a derivative with respect to the spacelike worldsheet co-ordinate. These condition will prove 
very useful in what follows. 
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Upon examining (|A8|I . we see that in order to make use of the spacetime delta functions, one would have to use 
the following transformation: 

So = dX ° dXa , (A13) 

\x°\\x ,a \ y ' 

where X a is the string co-ordinate field along any wound compact direction (which one we pick is insignificant). The 
term in the denominator comes from evaluating the jacobian of this transformation, subject to the condition ljA12f) . 
Note that we picked the particular co-ordinates X , because it is a monotonic function in time, and X a because it is 
easy to evaluate for wound strings. Using the constraints (|A10|) and IjAllfl . we see that: 



2na' / 1 
X° = -== , \^PiP 3 + 7 -—g ij X«X'i , (A14) 



V~9oo V {2ira'y 
where we use the fact that in the gauge we're in: 



P. = (A15) 

The expression in the sqaure root in (|A14|) is given by the Lq constraint in our constraint algebra [2j| (see also l|21|) 
and USD: an d is equal to the energy of the closed string: 



e = \ \Pd\ 2 + {n^n) + —(w :1 w) + — [2(n, w) + i(N - 1)], (A16) 
V a a 

where the worldsheet zero modes give us the contributions \pd\ for momentum along the non-compact directions, as 
well as the terms containing the winding and momentum quantum numbers for the compact dimensions. All the other 
fourier modes give us the oscillator contributions. We write the above as: 



A>° = 4^L e . (A17) 

V-300 

As for the second factor (|A13(I coming from the Jacobian, we see that for any string wound w times around the a th 
direction: 



X' a = w. (A18) 

Recall that we have to sum over all zeroes of the delta function along the a th direction when evaluating the integral 
(|A8|) after our change of variable (|A13|I . Hence the contribution 1|A17(I is cancelled by the fact that when a string 
winds w times around a particular direction, the argument of the delta function is zero precisely w times. Hence we 
evaluate (|A8p using the gauge fixing conditions and the results just obtained (|A14|) - (|A18|) as: 



q-,0 



rpl 

J i 

rp<t 



5 U -\X T -y T ), 



P Pi <rD-lf V T t\ 

-5 {X — y ), 



1 ( n l 



(A19) 
(A20) 
(A21) 



where g s is the determinant of the spatial part of the metric, and e defined by l|A16|) . We have deliberately only 
indicated the diagonal elements as they are the only elements which will survive a hydrodynamical averaging. We 
do this averaging as follows: keeping the quantum numbers p^,w,n and N fixed, we sum the contributions over a 
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distribution of such strings, where the momentum along the non-compact directions is distributed isotropically. Doing 
so gives us the result (I2*5l) - lj2"71) : 



P 



P 



Ho 

\fg~st' 



\Pd\ 2 /d, 



Mo ,n\ w 2 a b 2 



v^ 6 b l a ' 2 

This is exactly the same result one would obtain if one were to introduce the following action for the string gas: 



(A22) 
(A23) 

(A24) 



8 = - j d D+1 x^^^e . (A25) 

Note that in the above, we only assumed that the metric tensor was diagonal. Hence our result is valid even if we to 
begin in a frame such that goo was a function of time (rather than equal to —1), such as what would be required if 
we were to transform from the string frame to the Einstein frame using ()11JI . where we would like to end up with an 
Einstein frame metric of the form: 



9liV = diag[-l, a 2 {t),a 2 (t),a 2 {t),b\{t), ...b 2 p {t)]. (A26) 

Hence our energy momentum tensor in the Einstein frame would simply be the transformation of the above result 
using As for the issue of the energy momentum tensor of a D-string gas, it is straightforward to check that if 

we are in a regime where: 



°*(*°) « J_ (A27) 

then the only difference in the steps leading to the above derivation would be terms proportional to derivatives of 
<1> with respect to the fields X^. Since these will appear in conjunction with terms involving worldsheet derivatives 
of the fields A M , which will be of the string scale, the above condition will ensure that they will be negligable. For 
example if we demand that dx»& < lCP 6 /v / a', then almost certainly such terms will be insignificant and the net 
effect of having the extra dilaton coupling: 




(A28) 



will be the replacement of a' with a — e*/7i in (|A22|) - (|A24(I . This is essentially a generalization of the argument 
given in the appendix of [jj where the weak time dependence of the metric was shown to give a spectrum (|A16|I that 
evolved adiabatically (where the metric factors were allowed to depend on time). 



APPENDIX B: SUMMATION OVER MASSLESS MODES 

Consider the expression for the energy of a closed string l|24|l at the self dual radius, with the dilaton at its 
expectation value (j a b = a'5 a b, oj = 0): 



m 2 = 



1 r 



Q' 



{n + w,n + w) +A(N - 1) 



(Bl) 



We can read off that from this that the massless states are those which simultaneously satisfy the following set of 
equations: 



(n + w,n + w) = 4(1 - N), 
N +{n,w) > 0, 



(B2) 
(B3) 
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where the last equation is the level matching constraint. The quantum numbers that satisfy these conditions are 
tabulated below: 



N 


(n, n) 


(w, w) 


(n, w) 


1 











1 


1 


1 


-1 





1 


1 


1 





2 


2 








1 


3 








3 


1 








4 














4 






We note that whether or not the particular modes indicated above are realised or not depends on the number of 
extra dimensions. For instance for one or two extra dimensions, it will not be possible to satisfy (w, w) = 3 etc. 
We shall be most interested in the first four possibilities, classed separately above. This is because these states are 
not only massless at the self dual radius, with lj — 0, they are also massless to first order away from this point in 
moduli space. The first possibility (N — 1, n = w — 0) corresponds to unwound gravitons. As mentioned before, 
these states do not directly contribute to the dynamics of the moduli fields, as can be seen from l|41[l and (|42|) . The 
second possibility, in light of the level matching constraint can be seen to imply states of the form N = 1, n = ±u>, 
(w, w) = 1, corresponding to singly wound strings. These strings have an equal and opposite quantum of momentum 
along the same dimension, at oscillator level 1. The third possibility (again re-applying the level matching constraint) 
yields N = 0, n = w, (w, w) = 1 as massless quantum numbers corresponds to a singly wound state at oscillator level 
zero, with one quantum of momentum of the same sign along the same dimension. The fourth possibility (iV = 0, 
(n,n) = (id, w) = 2, (n,w) — 0) will also prove to be relevant. The second class of states are all at oscillator level 
zero, and correspond to various multiply wound/unwound strings with/without motion along various cycles of our 
torus. We now examine the behaviour of states with the quantum numbers tabulated above, as we move away from 
the rather special point in moduli space that we began in. To do this, we perturb our moduli as follows: 

7 =I-A, 
r l =I + A ; ||A||«1, 
oj = r) ; T) <C 1, 

where A is has a matrix norm that is much less than 1 . Expanding the formula for the mass of a closed string 124(1 , 
we see that: 

a'Sm 2 = (n,An) - (w,Aw) -4<q[(w,w) + (n,w) + 2(N - 1)], (B4) 

from which we see that only the first class outlined above remain massless to first order. For the first three states in 
the above table, the contribution from the metric fluctuations vanishes because ft and w are equal to each other up to 
a sign. For the states with quantum numbers given by N = 0, (n, n) = (w, w) = 2, (n, w) = 0, only those states that 
have entries in the same rows have vanishing fluctautions. That is, the condition (n, w) = implies that these states 
either have no entries in common, or two entries in common. The contribution from the metric fluctations would be 
non-vanishing in the former case, but for the latter, states of the form below remain massless to first order: 

w = (l,l)n = ±(l,-l), 
w = (1,-1) n = ±(1,1), 
w - (-1,1) n = ±(1,1), 
w - (-1,-1) n = ±(1,-1) 

where we only indicate the entries which the two vectors have in common. Our problem is now reduced to the following: 
having determined which states enter the summations in (|41() and (|42|l , we have to compute their contribution to the 
driving term taking care to account for all of them. We shall only compute the driving term for the radial moduli, as 
the case for the dilaton proceeds identically. Recalling that the unwound graviton states (N = l,n = w = 0) do not 



27 



contribute to the driving term, we compute the contribution of the states with quantum numbers N 
(w, w) = 1 to the driving term of (|41|l as: 



1, n = ±w, 



2 x 



8irG D fl _J_ 
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(B5) 



The factor of 2 comes from the overall degeneracy of the states which can appear either as, for example w — 
(0, 0, 0...1, 0...), n — (0,0,0... — 1,0...) or with the opposite signs, with the factor \pd\ coming from the factor of the 
energy in the denominator in the summand in (|41|) . The states with quantum numbers N = Q, n = w, (w,w) = 1 
sum to yield an identical contribution to the driving term. It is non-trivial, but straightforward to compute the 
contribution of the states corresponding to N = 0, (n, n) — (w, w) = 2, (n, w) = 0, where we find that we end up 
with the identical contribution as above, except now with the prefactor of 8(p — 1), where p — 1 has the combinatorial 
origin as being the number of ways one can pick two entries to be identical out of p choices, and 8 corresponds to the 
overall degeneracy of these states (as indicated above) . Hence the overall driving term is 
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(B6) 



where expanding around b a = VcV, u> = gives us (|43|l . The evaluation of the driving term for the dilaton proceeds 
analogously, as does the case for the summation over massless D-string modes. 
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